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Abstract. 

In order to study the origin of the strong Lya emission of high-redshift starburst galaxies we observed and modeled the emission 
of the z = 3.304 galaxy FDF-4691 (rest-frame EW]j,„ = 103 A). The observations show that FDF-4691 is a young starburst 
galaxy with a (for this redshift) typical metallicity. The broad, double-peaked profile of the Lycf emission line can be explained 
assuming a highly turbulent emission region in the inner part of the starburst galaxy, and a surrounding extended shell of low- 
density neutral gas with a normal dust/gas ratio and with Galactic dust properties. The detection of the Lyor emission line is 
explained by the intrinsic broad Lya emission and a low HI column density of the neutral shell. A low dust/gas ratio in the 
neutral shell is not needed to explain the strong Lya line. 
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' 1. Introduction 

I In all star-forming galaxies Lyff photons are produced by re- 

■ combination in HII regions ionized by young stars. However, 
these Lye photons are resonance scattered and thus have a large 
optical path in a neutral gas where they can be absorbed by dust 
grains (see e.g. Neufeld 1990). This may explain the absence 
of Lya emission from many local and medium-redshift star- 
burst galaxies (Chariot & Fall'1993'l. On the other hand, Lya 
emission is often the most conspicuous feature of the spectra 

\ of high-redshift galaxies. In fact, Lya emission is one of the 

■ most efficient tools to detect and identify high-redshift galaxies 
; (e.g. Rhoads et al. 200^ Hu et al.|l999j. Of 25 galaxies with 

z > 5 discussed by Taniguchi et al. (2003 1, all but two show 
strong Lya emission. Moreover, in short exposure spectra of 
high-redshift galaxies often no continuum is detected, while 
the Lya emission line is conspicuously present. Kudritzki et al. 
J2000i) argued that such galaxies must have a very low dust con- 
tent. However, as noted by Kunth et al. (1998) a higher escape 
probability of Lya photons could also be caused by a suitable 
velocity field, which reduces the number of resonance scatter- 
ing events. 

In order to clarify the cause of the strong Lya emission 
of high-redshift galaxies, we observed the Lya emission line 
galaxy FDF-4691 (from the catalog of Heidt et al. 2003 ) at low 
and medium spectral resolution. Moreover, we carried out ra- 
diative transfer model computations to reproduce the observed 

Send offprint requests to: C. Tapken, Heidelberg (e-mail: 
ctapken@lsw.uni-heidelberg.de) 

* Based on observations obtained at the ESO VLT at Cerro Paranal, 
Chile. 



complex Lya line profile and to constrain the velocity fields and 
physical conditions of the Lya emitting volume. Throughout 
this letter we adopt Qa = 0.7, Qm = 0.3 and Hq = 70 km 
s^'Mpc"'. 



2. Observations and data reduction 

The low-resolution spectrum {R = 200) of FDF-4691 (Figs. 
H] and E) has been taken from the atlas of Noll et al. J2004> . 
where the observational details are described. It is based on a 
total integration time of 550 min and has a continuum SNR of 
about 7 between 4000 and 8000 A. 

Medium-resolution spectroscopic observations of FDF- 
4691 were obtained in August 2002 in service mode, together 
with other objects. 8 exposures with 47 min integration time 
each were taken using FORS2 and the 1400V grism resulting 
in a total integration time of 6.25 h. The co-added spectrum 
covers the range from 4500 to 6100 A and has a resolution 
of R v 2000 (Av 150 km s '). The spectra were reduced 
using MIDAS. Since FDF-4691 is relatively faint (mi = 24.3 
mag), the continuum SNR of the medium resolution spectrum 
is < 3 throughout the observed spectral range. Hence this spec- 
trum provides little additional information on the absorption 
features. 

A much better SNR (< 60) was reached for the Lya emis- 
sion line profile. Fig.|3lshows the observed profile. The abscissa 
gives the radial velocity relative to the central absorption com- 
ponent (which corresponds to a redshift of z = 3.304). 
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Fig. 1. The low-resolution spectrum of FDF-4691. The dotted 
vertical lines indicate the positions of various expected lines. 
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Fig. 2. Scaled version of Fig. 1. The dashed line indicates the 
noise level. 



3. Results 

3.1. Basic observed properties 

Apart from the exceptionally strong Lya line the low -resolution 
spectrum (Fig. ^ shows a for this redshift normal starburst 
spectrum. The continuum break at the wavelength of Lye is 
caused by absorption of the Lya forest. The UV luminosity at 
1500 A is Luv = 1-63 + 0.06 x lO^"* W A"', which is typ- 
ical for galaxies between z = 3 and z = 4 in the FDF spec- 
troscopic survey (Noll et al. 2004). The rest-frame Lya equiv- 
alent width is 103 ± 15 A, while the line flux amounts to 
/Lya, = 18.8 + 0.6 X 10"^" W nr^. These are only lower Hmits 
since Lya photons are more affected by dust than the surround- 
ing continuum. The Lya luminosity is LLya- = 1.83 x 10-^^ W, if 
the emission is isotropic. 

As shown by Fig. |2l at least three prominent stellar wind 
lines are detected: NV ^^1239, 1243 and CIV ^^1548, 1551 
show P Cygni profiles and SiV AA1394, 1403 is mainly in ab- 
sorption. Hell ^1640 (total flux 5.3 + 1.6 x 10 W m ^) and 
CIII] ^1909 (total flux 8.8 ± 2.0 x lO^^i W m^^) appear in 
emission. A strong absorption {\EW\ = 3.8 + 0.9 A) blue-wards 
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Fig. 3. Lya line profile as derived from the medium-resolution 
spectrum of FDF-4691. The dashed line indicates the noise 
level. 

from the indicated position of CII /II 335 is discussed in Section 
3.3. 



3.2. AGN or starburst? 

Our Lya profile {FWZI ^ 2000 km s ') differs from Lya pro- 
files of high-redshift galaxies observed by, e.g., Dawson et 
al. (>2002j . which show narrower, asymmetric profiles with 
FWHM » 300 km s ' with a sharp blue cut-ofi" and a red wing. 
On the other hand, van Ojik et al. ( 1997 ) found for high-redshift 
radio galaxies symmetric Lya profiles with line widths of the 
order of 1500 km s As some of these profiles have tiie 
same line widths and similar profiles as FDF-4691, one might 
suspect that the Lya line of FDF-4691 is excited by an AGN. 

However, the hard radiation field of an AGN would be ex- 
pected to result in NV / Lya, CIV / Lya and CIII] /11909 / 
Lya line ratios which are not observed. For example, the ob- 
served CIV / Lya = 0.025 is ten times smaller than the ratio in 
the composite spectra of AGNs (CIV / Lya - 0.25, Osterbrock 
1989) . 

Moreover radio observations of the FDF have shown that at 
1.4 GHz no radio source > 100 yuJy is present at the position of 
FDF-4691 (Wagner 2003, private communication). 

Therefore, we conclude that the Lya emission of FDF-469 1 
is excited by the continuum radiation of the hot stars of this 
galaxy. 

3.3. Comparsion witli STARBURST99 models 

In order to constrain the starburst age and metallicity of FDF- 
469 1 we compared our spectrum with STARBURST99 models 
(Leitherer et al. 11999} . assuming a Calzetti et al. ilOOOi red- 



C. Tapken et al.: The nature of the Lya-emission region of FDF-4691 



3 




laOO 1300 1400 1500 

Wavelength [A] 



1600 



Fig. 4. Low-resolution spectrum of FDF-4691 (solid line) and 
STARBURST99 spectrum (dotted line). The discrepancy at the 
interstellar absorption lines (Sill, Ol/Sill and CII) is due to the 
fact that these are not modeled consistently with the extinction 
by STARBURST99. The dashed line indicates the noise level. 



dening curve. The best fit for the UV slope and the profiles of 
the stellar wind lines was found for continuous star-formation, 
low starburst age (< 20 Myr), LMC/SMC metallicity, a non- 
Salpeter IMF of » = 1.5 and E{B - V) = 0.1 mag (Fig.©. The 
metallicity derived from the CIV absorption equivalent width 
{\EW\ - 3.0 + 0.9 A) using the calibration given in Mehlert et 



al. J2002t is also in agreement with the LMC/SMC value. This 
metallicity is normal for this redshift (Mehlert et al. 2002). The 
reddening of E{B - V) - 0.1 mag is among the lowest ob- 
served in this redshift range in the FDF spectroscopic survey. 
Note that this value depends on the reddening curve used. The 
steeper extinction law of the SMC (Bouchet et al. 1985.) would 
yield E{B - V) - 0.03 mag. The low reddening is not unex- 
pected in view of the anticorrelation between the strength of 
the Lyo- line and the reddening as derived from the UV con- 
tinuum slope p found by Shapley et al. (2003 l. But deriving 
all the parameters mentioned above from UV data alone is not 
very constraining. Other combinations also provide good fits, 
including models with an instantaneous starburst and a very 
young age (^5 Myr). 

The CIV /1/11548, 1551 emission component is stronger in 
FDF-4691 than in the model. This may indicate the presence 
of early O stars or Wolf-Rayet stars, which are known to show 
such CIV profiles (Wu et al. .1983. Walborn & Panek .1984. 
Willis & Garmanv ll989> . The fact that the low-ionisation in- 
terstellar lines (Sill ^1260, Ol/Sill ^1303, CII /11335), visible 
in the STARBURST99 model, are not present in FDF-4691, 
is again consistent with the correlation found by Shapley et 
al. (2003) between the reddening and the strength of the inter- 
stellar absorption lines. The presence of the strong CII /11335 
line (blue-shifted by 1400 kms"' in respect to the Lyor Une) 
may be explained by the presence of WC stars in FDF-4691, 
which tend to show strong blue-shifted CII /11335 lines (Willis 
& Garmany[T989). WC stars may also provide an explanation 
for the strength of the Hell and CIII] emission. The presence of 
very young stars and WR stars together with the large equiva- 
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Fig. 5. Comparsion of the hya line of FDF-4691 (thin line) and 
the best-fit theoretical model (thick line). The intrinsic emis- 
sion profile is also shown (dotted line). 

lent width of Lyor indicate that the starburst is very young (< 5 
Myr). 

The star-formation rate of FDF-469 1 was obtained from the 
Lya luminosity (S FRi^ya = 16.5 + 0.5 Mg yr"') and from the 
UV luminosity (5 F/^uv - 17. 1 +0.8 Moyr"'). In both cases we 
used the calibration given in Kennicutt il998l l. The agreement 
indicates that Lya is not significantly attenuated relative to the 
continuum. 

4. Comparison with line profile models 

In order to constrain the physical properties of the Lya emitting 
volumes of FDF-469 1 we calculated model profiles using the 
finite element line formation code of Richling and Meinkohn 
(Richling et aL lMHI Meinkohn & Richling EUn^ . which is 
particularly well suited for calculating the radiative transfer in 
a non-static scattering medium. Since the width and the shape 
of the lower part of the observed Une profile is best explained 
by the emission of (high-velocity) turbulent H II regions, while 
the narrow central absorption can be produced only in a low- 
velocity H I layer in front of the emission region, we assumed 
a spherical two-component model with a central line emission 
region surrounded by a shell of neutral HI gas. Configurations 
of this type have been suggested and investigated earlier by 
Tenorio-Tagle et al. ( 1999 1, Ahn et al. (2001 2003 1 and Mas- 
Hesse et al. ( 2003 ). For the emission from the central region we 
assumed a broad Gaussian Lya emission profile. This assump- 
tion seems reasonable in view of the expected supernova rate 
and stellar wind activity in the observed compact starburst re- 
gion. The turbulent velocity of the central source was adjusted 
to fit the outer wings of the observed profile which are not af- 
fected by the low-velocity gas of the shell. 
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The computed model profile fitting the observed profile 
best is reproduced in Fig. 5. The theoretical model was con- 
volved with the instrumental profile. The parameters for this 
model are turbulence velocities of about 600 km s"' in the 
emission region, of 63 km s ' in the scattering shell, and an 
outflow velocity of the shell of 12 km s The central HI col- 
umn density of the shell was N(HI) w 4 x 10'^ cm"^, corre- 
sponding to an optical depth of the shell in the line center of 
TO ^ 5000. 

In the model described above the central absorption of the 
line profile is produced by the removal of Lya photons from the 
line center by multiple resonance scattering. Without dust ab- 
sorption the resonance scattering redistributes all these photons 
in velocity space to produce the blue and and red peaks fram- 
ing the central absorption. Calculations with Galactic dust/gas 
ratio and Galactic dust properties and without dust in the neu- 
tral shell showed that with the model parameters given above 
the dust in the shell had no detectable effect on the line pro- 
file. Dust in the central emission region, if not destroyed by the 
strong radiation field, only reduces the total emission without 
modifying the profile. 

In this model the strong Lya emission line is caused by an 
intrinsic broad Lya line, allowing a high fraction of Lya pho- 
tons to escape unaffected by the neutral shell, and a low column 
density of the neutral shell. Note that with Galactic dust/gas ra- 
tio the amount of dust in the shell is very low. 

Since the shell turbulence velocity appears too high with 
respect to the velocity of sound in (mainly) neutral interstellar 
matter and the rest intensity at the line center is somewhat to 
high in our model, we expect that we underestimated the neu- 
tral column density. A higher neutral column density would in- 
crease the separation between the two peaks, allowing to repro- 
duce a similar model as in Fig. [S] with lower shell turbulence ve- 
locity. Unfortunately, at present, our code cannot handle much 
higher optical depths. 

Additional model calculations with an increased dust/gas 
ratio and analytical models (Neufeld'1990'l suggest that a up to 
100 X higher column density of the neutral shell with Galactic 
dust/gas ratio will not result in a significant reduction of Lya 
photons in the shell. At higher neutral column densities the de- 
struction of Lya photons will become important. 

5. Conclusion 

Our observations have shown that the Lya galaxy FDF-4691 
is a young starburst galaxy with a (for this redshift) normal 
metallicity and a modest amount of reddening. The Lya flux of 
the galaxy appears not significantly more attenuated by dust 
absorption than the UV continuum. According to our mod- 
els the observed line profile can be explained assuming that a 
turbulence-produced broad profile of a central emission region 
is modified by frequency redistribution in a resonance scat- 
tering neutral shell around the central H II region. Although 
the model fit does not constrain the physical parameters of the 
model well (and the model may be much too simplistic), the 
computations demonstrate that the detection of the high Lya 
flux is caused by a low neutral column density and an intrinsic 



broad emission line. A low dust/gas ratio in the neutral shell is 
not needed to explain the strong Lya line. 
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